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THE STABH-ITY DlHIVAriVSS OF LOW-ASPECT-EATIO TRIANGULAR 
TONGS AT SUBSONIC AND SUPERSONIC SPEEDS 
By Herbert S. Hibner 


SUMMARY 


Low-arapect-ratio winss having triangular plan forme are treats 
on the assumption that the flow potentials in planes at right angle 
to the long axis of the airfoils are similar to the coiresponding 
two-dimensional potentials . Pressure distributions cauieed by down- 
Trard acceleration, pitching, rolling, sideslipping, and yawing are 
obtained for "vrlafTE with and without dihedx’al. The stability 
derivatives calctilated frean these distributions are expected to 
apply at both subsonic and supersonic speeds, with the exception 
of the •bransonlc region, up to a limiting speed at which the 
triangle is no longer narrow compared with 13ie Mach cone from its 
vertex . 


INTRODUCTION 


The aerodynamics of slender symmetrical pointed airfoils 
moving point foremost may be approximated as Munlc approximated the 
aerodynamics of slender airships (reference l) . For such bodies 
the flow is approximately two dimensional in plaiaos perpendicular 
to tlie axis of symmetry. The asounL’^ioii of two-dimensional flows 
leads to a very simple matlaematical procedure for obtaining the 
pressure distribution. Reference 2 Introduced this method and 
treated ihereby the slender pointed airfoil at an angle of attach:. 
The method is suited, as well, to the calculation of the pressure 
disti’ibutions due to normal acceleration, pitching, rolling, 
sideslipping, and yawing. In the present analysis the method in 
extended form is applied to the determination of these pressure 
distributions for a low-aspect-ratlo triangular plan form. The 
stability derivatives of the airfoiJ. are calcv.lated from these 
results . 
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The analyeis of the lifting airfoil shoved -ttiat if the airfoil 
is very slender (very lov aspect ratio) the restilts apply ■v/ell into 
the supersonic range i-rith no modification for the effect of 
corapressihility . The transonic region prohahly must he excluded. 

The stability derivatives of this report are expected to have a 
similar range of application. 

The principal part of this Investigation was carried out during 
March and April of 1946. 


SYMBOLS 


T 

x,yjZ 

u,v,w 






a 

3 

r 

AP 

p 

a 

h 

c 

c 

A 


flight velocity 

rectangular coordinates (fig. 2) 

incremental flight velocities along x~ ^ y- , 

and z-axes of figure 1, respectively; induced flow 
velocities along x-, y-, and z-axes of figure 2, 
respectively 

angular velocities about y-, and z-axes, 
respectively (fig. 1) 

component of velocity induced on upper surface 
parallel to stream velocity 

angle of attack 

angle of sideslip 

dihedral angle 

pressure difference between lover and upper surfaces 
of airfoil (positive in sense of lift) 

density of air 

seralwidth of triangle at distance x from vertex 

span (base of triangle) 

root chord (height of triangle) 

mean aero(3ynamio chord 

N 

aspect" ratio ; — ; 


5 s 


ib/2 


(local chord) ^dy 
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C 


G 

T1 

L 

L' 

M 

N 

Y 

F 



m 





edaa slope (| = ^ = f) 
area of triangle 

constant defined in equation ( 23 ) 
surface velocity potential 
cos"*^ y/a 
rolling moment 

normal force (approximately lift) 
pitching moment 
yawing moment 
lateral force 

suction force per tmit length, of edge 
lift coefficient 



rolling -moment coefficient 



pitching-moment coefficient 



yawing-noment coofficisnt 



lateral-force coefficient 


prof lie -drag coefficient 
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k vorticity 

Aq,Ai, . Fourier coefficients 

vjf induced surface velocity normal to wing leading edge 

s distance from -vrilng leading edge measured normal to 

edge 


Subscripts ; 

T.E. 

L.E. 

E 

L 


distance of center of gravity for^rard of ^o 


at trailing edge 
at leading edge 
at rigjit leading edge 
at left leading edge 
Subscripted parentheses; 

( contribution due to angle of attack 

( )p contribution due to dihedral 

I'Jhenever a, a, q, p, (3, and r £ire used as subscripts, a 

nondimensional derivative is indicated and this derivative is the 

slope throu{^ zero. For example, C n,^ 


"P 


time . 


^C. 






.Bp/ 


■^p=0 




BCjn 

; Cmq « 
d?=0 




'SOj' 




• 

r=0 


ntiation with respect to 


All angles are measured in radians 
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ANALYSIS 

SCOPE 


Tlie stability derivatives treated herein are listed, together 
■v^ith the values found, .for them, in table I. The derivations that 
follov give the values with reference to the principal "body axes 

of figure 1 -sfith origin at the aerodjTiamlc center ^^,0,0^. Con- 
version has been made to the system of stability axes shown in 

2 

figure 3 with origin a distance x^g ' ahead of iiie ^ point 

(transformation equations in reference 3) • Table I comprisos 
par all el columns ^rhich preseiit the values relative to both systems. 


GPSPAL 


Consider a slender isosceles triangle moving with its vertex 
foremost along its longitudinal axis, as in figure 1, with velocity V. 
Small linear disturbance velocities u, v, and w along the x-, y-, 
and z- axes and small angular disturbance velocities p, q, and r 
about the x-, y-, and z-axes, respectively, may bo contemplated. .Angle 
of attach gives rise to w, sideslip, to v^ and rolling, pitching, 
and yamng correspond to p, q, and r, x-espectlvely . 

As an example, suppose the sole disturbance velocity is w, 
caused by. .angle of attack. (Hiis case forms the subject of reference 2.) 
The triangular airfoil is assumed to be moving forward with velocity V 
and downward with the small velocity dV . The airfoil section is . 
assumed to be very thin; therefore, only the dovTnward motion disturbs 
the air. The triangular plan form is also assumed to be very slender 
so that the edges are nearly parallel. The flow in any. plane . 

X n constant (coordinate system of fig. 2) due, to the downward 
motion is .thus almost two dimensional. It may be expressed by the 
tvra- dimensional potential of a horizontal straight line moving, 
o-ownward with velocity dV . The horizontal straig^ht line is then 
the section of the airfoil cut by the plane x = constant. 

Planes x = constant may be taken anywhere from the apex to the 
trailing edge . , . . ' 
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In case the disturbance velocity is rate of roll p, the 
straight line is to be regarded as rotating with angular velocity p . 
The other cases are some^diat more complicated and are discussed in 
detail in subseguent- sections. In all cases ^ hoirever, the initial 
problem is the determination of the tiro- dimensional surface potential 
for the floi-r about a stral^t lj,ne with assigned boundary conditions . 

With the tvo- dimensional surface potential knoim for the motion 
of the section x = constant, the pressure difference between the 
upper and lower surfaces (positive upward) is obtained from 


AP = 2pVUq 


( 1 ) 


In this equation Uq is the component of the velocity induced on- the 
upper surface parallel to the; stream direction and is obtained by 
differentiation of the potential in the stream direction. Equation (l) 
expresses Bernoulli's law irith the approximation of small disturbances. 

The assumption that the triangular plan foim is very slender 
is expressed mathematically by the re3.ation C « 1. The qviantity C 
is the slope of the sides of the triangle relative to the stream 
direction and is equal to one-fourth the aspect ratio. The pressxrre 
distributions derived on this assumption can be shown to bo valid 
only to the first order in C. (See reference 4-, appendix- A.) Thus, 

2 A^ 

terms of order C neglected in comparieon with 

unity wherever they appear in the analysis. 

The validity of the analysis depends on the assumption that the 
disturbance parameters a, 3, pb/2T, gc/2T, and rb/2V are- 
small in comparison with unity. As in ordinary lifting-line irfng 
tlioory, however, terms in z? are of interest. Some such terms 
arise without approximation in the transformation from principal 
axes 'bo stability axes, but- others of order C (and hence not 
neg3.igible) result from retention of terms of order Cf^. For 
consistency, therefore, it has appeared necessary to, neglect all 
terms of order in the treatment. 

In the determination of certain of the stability derivatives, 
two cases will be considered. Case 1 is for a configuration having 
no dihedral and case 2, for a configuration having a small dihedral 
ang].e . 
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IDERIVATTVES Ci,^, Cj,^, AND 

The derivative is obtained ip. reference 2 . For 

accelerated motion the local pressure difference therein must be 
increased by the term 


= 2p II 

evaluated relative to axes fixed in the airfoil. This is 


AP = 2p l^a (2) 

da 

If the (email) angle of attack is a, ^e floir pattern in a 
plane cxxtting the airfoil at a distance • x frc«a the nose is the 
two-dimensional flow caused by a flat plate having the normal 
velocity aV . The surface potential is (reference 2) 


p » toJa sin ^ 


*= taV \ja^ - y^ (3) 

where cos t= 2, and the sign changes in going from the iipper 

GL 

surface to the lover surface of the airfoil. Differentiation of 0 
with respect to a, and substitution in equation (2) yield 


AP = 2paVa s3n 


Integration over the plan form gives the total incremental lift 
caused by it as 


L* = :^pVb^ca ‘ 
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This lift dlvlcLed "by incremental lift coefficient, 

and the derivative of this coefficient with respect to do/2V is 
the stability derivative C^. • It is 



( 4 ) 


The 

at 


center of pressvire of the distribution of AP 
X = ^ . The pitching moment about x = |c is 


is found to be 
therefore 



This moment divided by ^"V^Sc is the pitching“moment coefficient, 

and its derivative with respect to aic/2V is the stability 
derivative Cm . * It is 




(5) 


DERIVATIVES Ct and Cm 
■^q "“q 


An angular velocity of pitch q introduces a variation of 
angle of attack along the x-axis. 


a =* Oo 


' 1 °)' 


p 

where ccq is the angle of attack at x = ^.c . This variable a is 
to be substituted in equation (3) for the potential. 


0 = aV\j S - y^ 


= aVa sin ti 
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Eqimtion (l) for the pressure difference heti/een the upper and 
lower surfaces of the airfoil may be written 


AP = 2pV 

Carrying out the indicated differentiation gives 


AP = 2pT 


IToto + g 


(x-|c) 


CSC T ^ 2pVqa sin 


The integration of AP over the area of the triangle gives 
the value of lift found in reference 2 for an angle of attack Oq, 
pl\is the additional term. 


L' = § 

The coefficient is formed by division by —pV^S, and its derivative 
with respect to qc/2V is the stability- derivative *^Lg * 



(6) 


The integration of 


( 1 ° 


AP over the area of the triangle 


yields the pitching moment about the reference point 
moment is 




This 






2T 


The coefficient is defined as the moment divided by ip7®Sc and its 
derivative with respect to gc/2V is the stability derivative . 


The derivative is 
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DERIVATIVE Ci 
ip 

For the airfoil in rolling motion, any section x = constant 
is a rotating straight line. Lamb (reference 5) Gives the potential 
of the two-dimensional flow produced hy the rotation of an ellipse 
about its center. In the limiting case for which the ellipse 
becomes a straight line, the siirface potential is 


= |B7\la2 - (8) 

where p is the angular velocity, a is the semii’Tidth of the line, 

and cos ri = 2. . 

a 

Equation (l) for the pressure distribution takes the foim 

AP » 


By use of equation (8) , with C 


2pV 


2pV 


dx 

da dx 


da 

dx^ 


this expression becoraeB 


p 

^ =s pVpC X 



«= pVpC^x cot q 


(9) 


This antlsynmetric pressure distribution due to rolling was first 
obtained and shown graphically in reference 4. Figures 4 and 5 
herein are reproduced from this reference. Figiure 4 may be coEpared 
with figure .3 of reference 2 which shows the pressure distribution 
due to angle of attack. 
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The spanwlse loading, following reference 2, is merely 

^ = 2pY^6 

ay 


wl-fch 0 evaluated at the trailing edge of the airfoil. The die- 
trihution is sho^m. in figure 5 « 


The integration of 
moment 


dl.* 

ay 


■y ay 


across the span gives the rolling 


I. 



Division hy ^V^Sh converts this moment to coefficient form, 

and the derivative with respect to ph/2V ' is the stability 
derivative Ci . It is 


C 


‘'P 


3t b^ 

32 S ■ 



* ( 10 ) 


DEBIVAnVE C 7 
Case 1, No Dihedral 


The pressure difference across a thin airfoil in steady flow 
has been given by equation (l) 


AP «= 2pTUq 


where is the con 5 >onen,t of induced velocity parallel to the 

stream direction, measured along the upper surface. If sideslip 
occtirs, the stream direction is inclined relative to the x-azis of 
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the airfoil hy the sideslip angle p. Thus, with p positive for 
sideslip In the positive direction along the y-axis, 


Uq e cos P - sin p 
^ ox dy 


<« 




where 


p^« 1. 


Hence, 


a? =2pv(g- (11) 

The surface potential 0 for iiie diaturhance velocity depends 
only on the normal velocity of points of the surface. The potential 
is therefore unaffected hy sideslip when there is no dihedral. In 
the present case the normal velocity is oV due to an£^e of attack 
and the appropriate potential is that discussed in the section 
on and Cjjj^ and given by eqiiation (3) • Carrying out the 

differentiations indicated in equation (ll) gives, after 
simplificaticaa. 


« 2p'V^a CSC T] + p cot 


( 12 ) 


The symmetric first term is the lift distribution in the 
absence of sideslip. Hiis term yields no rolling moment. The 
antisymmetric second term contributes the rolling moment 


The coefficient is formed by division by ^V^Sb, and its 
derivative with respect to p is the stability derivative 
It is 



( 13 ) 
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Case 2, Dihedral Angle F 

If the tring has a email dihecLral angle P (case 2), the ang3.e 
of attack on the left panel Is reduced in sideslip "by the amount PP 
and the angle of attack on the ri^t panel is increased hy that 
amount. The flow pattern in a plane cutting the airfoil at a 
distance x. from i3ie nose can he obtained hy a slight modification 
of tlie classical thin airfoil theory. (See reference 6.) The 
left edge and rlgjht edge of the section at x are to he identified, 
respectirely, with the loading edge and trailing edge of the section 
of thin airfoil theory. !Ehe section is regarded as a small 
deviation from its chord. A diatrihutlon of vorticity along the 
chord of the section is imagined. Paraphrasing Glauert (reference 6), 
the induced velocity w is determined for points on ilie chord hut 
may he taken to ha the same for the corresponding points of the 
section Itself . The direction of the resultant velocity adjacent 
to the airfoil must he parallel to the surface so that at each 
point of the left panel 


2 = -|3r + a (l4a) 

and at each point of the right panel 

I = PP + a (l4h) 


The potential corresponding to a is already known; therefore, only 
the case for cc = 0 need he treated. 

The vorticity assumed in the thin airfoil theory (reference 6) 
has a net circulation to satisfy the Kutta condition at the trailing 

edge. The addition of a term -2 t(aq + ^ 1^ esc T) (with 17 written 

for Glauert 's 0) eliminates the circtilation while retaining certain 
mathematical properties. With this added term the vorticily is 
given as 


k 


2V(Aq cot T] 


■^1 

CSC n 



(15) 


1 
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Tho velocity potential on the upper surface of the section is 
related to the vortlcity k hy 



The angle t], originally Identified with Glauert's B, is now 
defined differently so that y = a cos . Then Idle integration 
of the equation for 0 with dy *= -a sin t) dq gives 


0 a av{- 


Aq sin 11 + sin 2q 


00 


& 

2 


sin fn + l)n sin (n l)n 1 

n+1 " n-1 ' 


Equation (l6) expresses the upper surface potential for an 
arhitrary distribution of induced vertical velocity along a line 
in two-dimensional flow irlthout circulation. 

The coefficients in equation (l6) are still to be evaluated. 
The calculation given on pages 88 to 90 of reference 6 when 
applied to equation (15) leads to 


^ = “Aq + 2_ cos mi 
1 

for the ratio of the induced downward velocity to the stream 
velocity. The coefficients are given by tho theory of Pourler 
series as 



(IT) 
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For the sideBlipping airfoil section with dihedral angle, the 
■boundary conditions on w/V has "been given in ©qiiations (l4a) 
and (ite) . By eqtiation (l?) the coefficients in this case are, 
for a •= 0, 


Aq « 0 


An 


^ Ag£, 

Tin 


. nit 
sin -g- 




(18) 


With tlie coefficients given in equation (l3), then, equation (l6) 
represents the additional potential duo to dihedral that may ho 

substituted in equation (U) * The terra p is found to he of 
2 

the order 3 and thuis may Ije neglected for the present purpose. 
Integrating tlie pressure chordwise gives the incremental spanwlse 
load dlstrihution caused hy dihedral as _ 
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By equation (l8) A]_ 

edge . a o Thus, 


s and Ao = - , and at the trailing 

jt ^ 3rt 

for the present case 




The coefficient is formed hy division by ^v^Sb and its derivative 
■vhLth respect to j3 is the stability derivative 0^^ for zero angle 
of attack. The derivative is 


■■ 

The con^lete stability derivative is obtained by adding to tho 
preceding equation tho contribution of angD.e of attack given in 
equation (13) so that 



(19) 


DERIVATIVE C? 

‘t 

Case 1, Ho Dihedral 


Tlae reference point for these calctilations is at, a distance ^ 

from the vertex of the triangle, measured along the x-axis. Let the 
stream velocity at this point be V. Then, if the yawing velocity 
is r, the longitudinal velocity at (x,y) is V - ry and the 


sideslip velocity is -r 


- ¥)■ 


Let these expressions replace V 


and |3V, respectively, in equation (ll) so that 


(V - ry) 



Z!$> = 2p 
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The appropriate potential 0 i.g that given "by eqtiatlon (3) , 
vhich is 

0 *= dVa sin Tj 


Carrying out the indicated operation gives 


^ = 2p 


V^cxC cac I) - rVa^ “ ^ 


cot TJ 


( 20 ) 


whore C = ^ « 1. The term aC = aoC^ end is negligible in 
dn 

comparison with x to the firat order in C» 

The symmetric cooecant term is the lift distribution in 
straight fli^^t; it yields no rol3.1ng moment. The antisymmetric 
cotangent term contributes the rolling mamont 


Tc na 

L a 2p I rVa 
Jo J-a 


(. - |c) 


cot TJ 7 dx dy 


With T] = arc cos this moment is eTOluoted as 


L = -—pvretb^c^ 
144 


i^V^Sb, 


and its 


The coefficient is defined as the moment divided by 

71 

derivative with reeuect to rb/2v is the stability derivative C-, 

h 

for zero dihedral. It is . - ^ 
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Case 2^ Dihedral Angle r 


In the section on Ci 


dihedral was shown to add an additional 


term, which is given hy equation (l6), to the potential, where the 
constants are given hy equation (l8) . In the present section 

-f(^ " 3 ?) 


the variable sideslip angle 
sideslip angle 3 of the section on 
to the additional term is thus 


renlaoes the constant 


The rolling moment duo 


|pv2a3(Ai - Aj) 


T.S. 


as was found in the section on equation (l8) , 

with P = - — /x - 

V\ 3 / 




The substitution of -the values of a, Aj_, and A^ at the trailing 
edge, where = g" x = c, gives 

L = i^Trfb3c 

The coefficient is formed as before by division by ^pT^b, and its 
derivative with respect to rb/2V is 



2 

9 


r 
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Hiis result must "be added to the result for case 1 (no diheiSral) 
to give the total value of the staMlity derivative, as follows: 


C, =* SS + ■ ^ 22 ) 

Ir 9A 9 


DiffilVATmilS Cy 


AND 



Case 1, ITo Dihedral 


The side force and jawing jncssent relative to tody axes are 
contrihutod entirely hy suction along the leading edge of the wing. 
This suction may he evaluated hy considering the triangular wing to 
have a small thlokness so that the sections x = constant are 
ellipses. Tlie lateral component of the pressure dlstrihution is 
determined and integrated, ^io aTc^proach is given in detail in 
reference 4 for the case of rolling at zero aiigle of attack, and 
its extension to the case of i-olling with a small anf^Le of attack 
was made in the original determination of Oy^ and in the 

present analysis. A very much simpler method of evaluating ttie 
suction is s^lggested in reference 7> howover, and this method is 
adopted herein. 

Consider a condition for which the induced surface velocity 
normal to the edge is of the foim 


Ylf = (23) 


in the iimiediate noi^horhood of tlxe edge, XThere s is distance ^ 
from the edge and G- is a constant. Reference 7 points out that 
for such a flow there is a suction force per unit length of edge 
which is 



(24) 


in an incompressihle fluid, ^or the narrow triangles discussed in 
this paper the component of the stream velocity normal to the edges 
is inherently small corroarod \rLth the velocity of sound over the 
range of stream Mach numbers considered. TIuis, no ccmpressihillty 
correction is necessary. 
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For the triangular TriLng in rolling motion the Induced velocity 
component u has been obtained as a factor in eqtiation (9) end may 
he -written 


^1 



Angle of attack -?rill give the additional contrihutlon (reference 2) 



The total induced, velocity component u on the upper surface is thus 





Very near -the edge -this escprossion is approximately 



Trhere -fche plus sign refers -to -the ri^t edge and -the minus sign -to 
•the left edge. 


If a similar calculation is made for v 


it is found -iliat 

.»jf 

velocity Vu^ + v^ 


as -the edge is approached -the resultant induced 
hecomes normal to -the edgsi Thus, -the normal velocity near the 
edge is 
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Tjj = 


s/ 


1 + 0 “= 


■u 



to the first order In C. The perpendicular distance of point (x,y) 
from the edge is, to the same degree of accuracy. 


= III} 

The induced surface velocity very near the edge my therefore he 
expressed approximately as 





•vdilch Is of the form of equation ( 23 ) already considered. The 
suction force per tmit length of edge by equation (24) is thxis 

F B «pG-^ 

= ^Cx:^Y® + ^^C^x? ± aVpC^ 


where the plus sign refers to the ri^t edge and the minus sign 
refers to the loft edge. 

The lateral component of this suction force is given by 




no' 

y = J - Fl) dx 

= 2pC^c3aVp 
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The coefficient is formed hy dlTislon hy 
with respect to ph/2V is the stability 


3:o'V^S, and the 
derivative Cy . 


derivative 
It is 


(%)a “ 


= ~jta 


( 25 ) 


The moment of the leading?, -edge suction about tha vortex of the 
triangle is approximately, for c2 « 1, 


' Jo ^ 


= -^C^o^dVp 


The moment about the reference point 


(|c,o) 


is 


N » Wg + |cY 


= -■^pC^c^aVp 


The moment coefficient is formed by division by ^V^Sb, and the 
derivative with respect to pb/2V is Idle stability derivative C 


"P* 


It is 



( 26 ) 


Case 2, Dihedral Angle F 

To the first order in the dihedral angle F , dihedral will not 
change the pressure distribution. The inclination of -the wing panels, 
however, will give rise to a lateral force component, as follows: 
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y 




^ dx 


cly 


for the right Tring panel, and a similar expression vilii opposite 
sign for the left wing panel. The presstire difference ^ has 
heen evaluated in equation (9) of the section on With this 

value the integration gives 


Y = -^vrb^p 


Division hy 

the increment to the 


and differentiation with respect to ph/2V give 
stahillty derivative Oy catised by dihedral as 



By addition of the value obtained for ihe case of no dihedral, the 
complete derivative Oy is 



Irto- ■ 

3 3 


( 2 ?) 


The pressure distribution is such that along any radial line 
from the vertex the pressure increases in proportion to x. For 

such a distribution the center of pressure on each panel is at ^ 

/e 

from the vertex, 
is 


The yawing moment about the reference point 
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Division "by ^V^Sh and differentiation with respect to ph/2V 
give the increment to caused hy dihedral as 

(%)r “ 

By addition of the value obtained for no dihedral idle ccmipleto 
derivative is therefore _ . 


S ” “9A 

DERIVATIVES Cy AND 
Case i. No Dihedral 


A little sideslip can readily he shown to have negligible 
effect on ihe symmetrical distribution of suction along most otthe 
leading edge . Near the trailing edge some modification may bo 
expected at subsonic speeds due to the altered direction of the 
trailing vortex sheet. Any lateral force and yawing moment would 
have to come from the small, disturbed region. Examination indicates 
that such a force or moment would bo of order o^p and hence zero 
to the first order in a. 


Case 2^ Dihedral Angle r 

The contribution of dihedral to the velocity conponont u 
induced in sideslip is obtained by differentiating equation (l6) 
with respect to x. A term dVC esc ti must be added for the 
effect of angle of attack. After insertion of the constants from 
equation (l8) the total velocity component u is obtained as 


u = oVC CSC T| + — prCV-s + cos 2tj cot tj 




rut 

2 


sin (n + l)n + sin (n - l)t ^ 
n + 1 n - 1 J 


( 29 ) 
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The only terms that approach Infinity near the edges (that is, 
as T]->0 or «) are the cosecant term and the cotangent term. In 

that region both terms hehave like appropriate sign, 

vhere s is the perpendiciolar distance from the edge. The velocity 
component u there may thus he -written 



where the plus sign refers -bo the ri^t edge and -the minixs sign to 
•the left edge. 

If a similar calcialatlon is made for v = it is found, as 

oy 

before that near the edge -the normal induced velocity Vj^ is 
approximately equal -bo u/C to "the first order in C. Th-us, vjj is 

of "the form and the corres-Donding suction force per unit leng-th 

edge is F «=■ as indicated in a preceding section. Substi-butlon 

of -the expression for G-, neglect of -terms of the second order in p, 
and simplification gives 

F » |pV^Cx(^ t (30) 


The lateral component of -this suction force is 

Yi « J* (Fr - Fp,) dx 

= 2p-V®Cc^apr (31) 


There is an additional lateral force due -bo a component of 
•the pressure acting on the inclined panels. The part of the 
pressure caused by dihedral -vrlll con-tribute -terms of order 
To the first order in then, only -the pressure distribution 
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in the absence of dihedral (equation (l2)) need he considered. 
Further, only the antisymmetric cotangent term will contribute to 
the lateral force. The incremental lateral force is thus 


Yg - -2r / r ! ki|pv2|3 


JJ 


cot iq dy dx 


no 


= -4pT2apr 

UO U 


n«/2 


Ox cos T] dT) dx 


-2pV^Cc^apr 


(32) 


The total lateral force to the first order in F is + Y2 • 
This sum is seen to be zero, and hence the lateral force 
derivative Gyp is likewj.se zero. 


Equation (12) shows that sideslip gives rise to a pressure dis- 
tribution that is constant along radial lines frcaa the vertex of the 
triangle . (Such a pressiare distribution defines a conical f lovr 
field.) Equation (12) is for- zero dihedral, but equation (16) leads 
to the same behavior for the triangle with dihedral. The center of 

pressure on each panel of the -briangle will be on the line x 
So also will the center of pressure of iiio leading-edge suction. 
There is thus no yawing moment about the reference point 

the stability derivative is zero. 


2 

3° 


and 


DERIVATIVES Cy^ A3SD 


In the case of yawing motion the local sideslip velocity 
is ~ suction F per unit length of the leading edge 

is obtained by substituting this local sideslip velocity for pv 
in equation (30) derived- for sideslip, as follows: 
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F = |pV^Cx 

(TThen p Tarles with x, as in the present case, equation (29) 

requires an additional term in proportional to F , This term 

is finite at the edges and, therefore, does not contribute to the 
expression for Y.) 

The lateral component of this suction force is 

Yl.“ J - I’D 




(33) 


The antisyme+ric cotangent part of the pressure distribution 
in equation (20) may contribute a lateral force because of the 
inclination of the panels . To the first order in F, this is the 
only contribution, according to the reasoning in the section 
on Oy and ^np* contribution is, to the first order in C, 


ij^jVorF j ~ 

Jo Jo 


dx 


4pVarF j j V " 37 


Cx coo rj dq dx 


10 ».'o 


(34) 


The total lateral force caiised by yawing 3Botion is + Y2* 
The lateral force derivative Cy is accordingly zero. 

The leading-edge suction gives rise to a pure yawing couple . 
This couple is conveniently obtained by conq^uting the moment of 
the suction force about the vertex of the triangle, as follows 
(C^ « 1) t 
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Jo 

= ^v^cocrcc^ (35) 


The pressure distrlhution on the inclined panels of the wing 
will have a lateral component that likewise contributes a pure 
yawing couple. Fran equation ( 2 o) the couple is^ to the first 
order in C, 


pc pa 

N 2 » -4rpVa II fx~ |c j cot T) X dy dx 

Jo Jo ‘ 

pc prt/2 

= “4rpToC J J ~ T) dq dx 


sx 


-^VorrCc^ 


(36) 


A third yawing couple irtll bo contributed by slcJn friction. 
(Skin friction has not been considered in evaluating the other 
stability derivatives because its direct contribution is expected 
to bo unimportant. ■■ The indirect effect of skin friction in 
Influencing idle pressure distribution via -ttio boundary layer may 
indeed be Important, although it is not treated.) 

The skin friction couple is approximately given by: 


Wo 





P + y 


dx dy 


where Vjj is -Hie resxiltant velocity and ■vp = (y 


3 is the local sideslip angle and equals 


E 


ry)^ H- 


•(x - 1 =) 


r^^x " 


and C 


-y -Do 

the section drag coefficient which is taken equal to wing profile 
drag coefficient. 


is 
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To the first order in r this expression is 


^3 “ 


nc oa 

f 

-w 0 


1 f " 

Vr 

(x “ |c) + 2 y 2 

- V^y " 

Jo J-a 


\ 3 y 



•^VrCD^Cc^d + 6c2) 


The total yawing moment caused hy yawing motion is Nj + N 2 ^ 3 ^ 
which is Just N^ hecause N^_ and N 2 cancel. Ihe coefficient is N-^ 

divided hy |o'V^Sb and the derivative with respect to ~ is tlie 
stahility derivative Cn^, . Carrying out these operations gives 



(37) 


A similar calculation shoTO that the side force due to skin 
friction is zero. 


EESULTS AND DISCUSSION 


The values obtained for the stahility derivatives are 
summarized in table I with respect to two systems of axes. One 
system is the principal body axes of figuK*© 1 with origin at the 


aerodynamic center 


• (§c,o,o) 


The other system is the stability 


axes shown in figure 3 tritla origin a distance ahead of 


the |c point. 


These stability derivatives apply to an isolated triangular 
wing in the limiting case of aspect ratio approaching aero. Applica- 
bility decreases \ 7 ith increasing aspect ratio, and an aspect ratio 0.5 
is estimated as the upper limit of utilitj'. The mathematical validity 
at the very low aspect ratios may be offset, perhaps, by error duo to 
the neglected boundary layer. 
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The argimentB for the effects of compressihility presented in 
reference 2 can he carried over to the present "srork. The stahility 
derivatives presented herein, therefore, are expected to apply at 
"both subsonic and supersonic speeds, with the exception of the 
transonic region, up to a limiting speed at which the triangle is 
no longer narrow compared with the Mach cone from its vertex. 

The over- all pitching-moment derivatives should he little 
affected hy the addition of a fuselage, the nose of ■vdiich does not 
project much heyond the vertex of the triangular \r±ng* The wing 
will orient the flow along the axis of the fuselage and thereby 
■vriLll eliminate much of- the unstable pitching moment of the fuselage . 
The flow will continue to be essentially axial along the part of 
the fuselage behind the wing because the low aspect ratio yields a 
doiTn-sTash angle substantially equal to the angle of attack. 

Theoretical- considerations suggest that the unstable yawing 

moment of the fuselage will add to the value of for Ihe wing 

p 

alone. Little effect on or is expected. Little effect 

on the rolling-moment derivatives is expected if the wing is mounted 
centrally on tlie fuselage. Hig^-wing or low-wing arrangements, 
however, should have pronounced effects on the effective dihedral. 
These conclusions are only tentative, and a proper evaluation of 
the ■^'d.ng-fxiseleige interference mtist be the subject of farther 
investigation. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., July 15, 1947 
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Figure 2,- Axes and notation used in analysis. 



Figure 3»- Velocities, forces, and moments relative to 
stability axes vrith origin at |c x^g. Principal 

axes of figure 1 dotted in for comparison* 
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Figtire it»- Distribution of pressure difference caused by rolling. 





Figure Distribution along span of normal force caused by rolling. 


